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We have carried out bound-state and low-energy quantum scattering calculations on He -f NH 
(^E~) in magnetic fields, with the NH molecule in its n = 1 rotationally excited states. We have 
explored the pattern of levels as a function of magnetic field and identified the nearly good quantum 
numbers in different regimes. We have used the bound-state calculations to locate low-energy 
Feshbach resonances. When the magnetic field is used to tune across such a resonance, the real 
and imaginary part of the scattering length show asymmetric oscillations and peaks with amplitude 
between 1 and 3 A. The scattering length does not pass through a pole at resonance. The resonant 
behavior is characterized by a complex resonant scattering length Ores- The corresponding inelastic 
cross sections show troughs as well as peaks near resonance. This may be important for efforts to 
achieve evaporative and sympathetic cooling for molecules, because it offers the hope that inelastic 
trap losses can be reduced by tuning close to a Feshbach resonance. 

PACS numbers: 34.50.-s,34.10.-|-x,03.65.Nk,82.20.Xr,34.30.+h 
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I. INTRODUCTION 

Collisions between ultracold atoms can be controlled 
by tuning the atomic interactions using applied magnetic 
fields Such techniques have been used to cause con- 

trolled implosion of Bose-Einstein condensates 0] and to 
produce dimers of both bosonic 0, H, 0] and fermionic 
[1, [^, [13, mi alkali metal atoms. Long-lived molecular 
Bose-Einstein condensates of fermion dimers have been 
produced [H, [ll [11, and the first signatures of ultra- 
cold triatomic and tetraatomic [1^1 molecules have 
been observed. It is proving possible to move molecules 
selectively between quantum states by either fast or slow 
sweeps of magnetic fields across avoided crossings be- 
tween bound states [Il| ■ 

The overall strength of the interaction between a pai r 
of atoms is characterized by the scattering length a [l9| . 
An important feature of elastic scattering in ultracold 
atomic gases is that a passes through a pole as the mag- 
netic field B is swept across a Feshbach resonance at 
constant kinetic energy [20[, 



a{B) = abg 



B — B„ 



(1) 



where Obg is a near-constant background scattering 
length and Sres and Ab are the position and width of 
the resonance. The scattering length can thus be tuned 
to any desired value, positive or negative. Positive values 
correspond to interactions that are overall repulsive and 
negative values to interactions that are overall attractive. 
The elastic cross section is given by 



1 + k^a^ ' 



(2) 



where the kinetic energy is i?kin = ?i^fc^/(2/i) and ^ is the 
reduced mass for the collision. The elastic cross section 
thus passes through a peak of height 47r/fc^ at resonance. 

It has recently become possible to cool molecules di- 
rectly from high temperature to the millikelvin regime, 
using methods such as buffer-gas cooling [2li, [^ 
and Stark deceleration [2^ . [25j . Polar molecules such as 
ND3 and OH have been successfully trapped at temper- 
atures around 10 mK, and there are a variety of propos- 
als for ways to cool them further, including evaporative 
coo ling , sy mpathetic cooling and cavity-assisted cooling 
[26I . [27l . [28| . Very recently, NH has been trapped at tem- 
peratures around 0.7 K by buffer-gas cooling in cryogenic 
helium. 

In previous work [23| , we have explored the possibility 
of controlling molecular interactions in the same way as 
atomic interactions. We have generalised the BOUND 
[13] and MOLSCAT packages to carry out bound- 
state calculations and quantum scattering calculations in 
applied magnetic fields for systems made up of a ■^■^+^1] 
molecule and a structureless atom. 

In our initial calculations on He + NH ('^E) [l^l, we 
used BOUND to locate magnetic fields at which bound 
states cross open-channel thresholds. We then used 
MOLSCAT to characterize the resulting low-energy Fes- 
hbach resonances as a function of magnetic field. For a 
resonance at which a bound state crossed the lowest open- 
channel threshold, we observed a pole in the scattering 
length that followed Eq. [TJ However, for a resonance in 
which a state crossed a higher threshold, we observed 
quite different behavior. The scattering length showed 
only a weak oscillation instead of a pole. The suppres- 
sion of the pole was attributed to inelastic effects. The 
calculations were for NH in its lowest rotational state 
(n = 0), for which inelastic coupling is very weak, and 
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the resonances were very narrow, but even so the amph- 
tude of the oscillation in a{B) was only about 9 A. 

A full derivation of the resonant behavior of the scat- 
tering length in the presence of inelastic effects has been 
given previously (32i | , so only a brief version will be given 
here to explain the basic physics and define notation. 
In the presence of inelastic collisions, the scattering ma- 
trix S that describes the collision in quantum-mechanical 
terms has elements Sn' . The diagonal S-matrix element 
in the incoming channel has magnitude |S'oo| < 1 and 
may be written in terms of a complex phase shift Sq with 
a positive imaginary part [33| . 



Sooiko) 



(3) 



where /cq is the wave vector in the incoming channel. This 
can be expressed in terms of a complex energy-dependent 
scattering length, a{ko) = a(fco) — i/3(fco) [M,!!^!, defined 
as 



^^^^^ _ -tan(5o(fco) _ J_ / 1 - 5'oo(fco) 
fco iko \l + Soo{ko) 



(4) 



a(fco) becomes constant at limitingly low kinetic energy. 
The elastic and total inelastic cross sections are exactly 

M 



and 



crol(fco) 



47r|a|' 



l + A:2|a|2 + 2fco/3 



ko(l 



2kol3) ■ 



(5) 



(6) 



A scattering resonance is most simply characterized in 
terms of the S-matrix eigenphase sum E, which is the 
sum of phase shifts obtained from the eigenvalues of the 
S matrix [s^, [sll . The eigenphase sum is a real quantity, 
and across a resonance follows the Breit-Wigner form, 



E(i?) 



tan 



2(-Brcs — B) 



(7) 



where Ebg is a slowly varying background term, B^^s 
is the resonance position and Tb is a resonance width 
(not the same as As in Eq. [J). The B subscripts in- 
dicate that we are considering the resonance as a func- 
tion of magnetic field rather than energy. The individual 
S-matrix elements describe circles in the complex plane 



Sii'{B) = S'bg, 



igBi9Bi 



B - B,,, + irB/2' 



(8) 



where gsi is complex. The radius of the circle in Sa is 
ISsi/^sl- T^^'^'^ partial width TBi for channel i is usually 
defined as a real quantity, but here we also need a corre- 
sponding phase (jji to describe the orientation of the circle 
in the complex plane, = Fsje^"^'. The width Tb and 
partial widths Tbi are signed quantities, positive if the 
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FIG. 1: (Color online). The small circles in the elastic S- 
matrix elements in the presence of inelastic scattering for He 
-I- NH (n — 0) scattering at iSkin = 10~^ K (green, smaller 
circle) and 4 x 10~® K (red, larger circle) Reproduced with 
permission from ref. [29l . 



resonant state tunes downwards across the threshold as 
a function of B and negative if it tunes upwards. For a 
narrow resonance, the total width is just the sum of the 
partial widths. 



Bi- 



(9) 



The partial widths for elastic channels (degenerate 
with the incoming channel) are proportional to fcg at low 
energy. We may define a reduced partial width jbo for 
the incoming channel by 



T^Boiko) 



(10) 



and the reduced width is independent of fco at low energy 
(typically below i?kin = 1 mK). By contrast, the par- 
tial widths for inelastic channels depend on open-channel 
wavefunctions with large wave vectors ki and are effec- 
tively independent of fcg in the ultracold regime. If the 
inelastic partial widths r^i are non-zero, they eventu- 
ally dominate r^o as fco decreases. The radius of the 
circle ([5]) described by S'oo thus drops linearly to zero as 
fco decreases, as shown in Figure [1] This is qualitatively 
different from the behavior in the absence of inelastic 
channels, where the circle has radius 1 even at limitingly 
low energy. 

The radius of the circle in S'oo is Tbo/Tb- For small fco, 
where Eq.fTOl applies, this is approximately 2fco7Bo/r'B'''- 
The formula followed by the complex scattering length is 



a{B) = flbg 



2iB-B,,,)/r^^' + i 



(11) 



where F'J°' is the energy-independent part of F^ (omit- 
ting Tbq) and ares is a resonant scattering length that 
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Magnetic Field (G) 



FIG. 2: (Color online). Real and imaginary parts of the 
scattering length across a Feshbach resonance in He + NH 
(n = 0), showing a small symmetrical oscillation in the real 
part (red) and a peak in the imaginary part (green). Repro- 
duced with permission from ref. |32| . 
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characterizes the strength of the resonance, 

_ 27-BO 2i(0o+fcoabg) 
"res — ijj^i e 

^ B 



(12) 



Both Arcs and the background term Obg can in general be 
complex and are independent of kg at low energy. How- 
ever, in the special case where the background scattering 
is purely elastic (obg is real), unitarity requires that the 
circle in 6*00 must loop towards the origin as shown in the 
upper panel of Fig. [31 This requires that a^cs is also real. 
Across the width of the resonance, the real part a{B) of 
the scattering length a{B) then oscillates about Cbg by 
±ares/2 and the imaginary part peaks at 0(B) = arcs- 
This was the behavior seen for a{B) in ref. "2!? for He + 
NH('^S) collisions with NH in n = states, shown in Fig. 

m 

The purpose of the present work is to explore a rather 
more complicated case, with significant background in- 
elastic scattering. Eq.[TT]still holds, but arcs can be com- 
plex and the circle in 5*00 then does not point directly 
towards the origin. This behavior is shown in the lower 
panel of Fig. [31 The elastic cross section is given by 



^oi(s) = -j|i-5oo(fco)r 



(13) 



so at any value of B it depends on the distance between 
Sqo and the point X at 6*00 = +1. However, the total 
inelastic cross section is given by 



^inel 



^(l-l^oo(MP) 



(14) 



and thus depends on the distance of 6*00 from the unit 
circle. If the circle in Sqq does not point directly towards 
the origin, it is clear from Fig. [31 that the total inelastic 
cross section can show a trough as well as a peak near 
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FIG. 3: The distinction between an S-matrix circle in the 
complex plane that points directly towards the origin O (top) 
and one that does not (bottom), showing the relationship be- 
tween dres and cj)o. 



resonance. This offers the hope that resonances can be 
used to reduce inelastic rates as well as increase them. 

In the general case, the explicit exp ressions for the real 
and imaginary parts of a{B) are |32| 



a{B) = abg + 
P(B) = /3bg + 



arcs [2(5 - Bres)/rg°'] - A 

5rcs)/r*iJi' + l 
s [2(5 - Brcs)/r^ 



[2(5- 



incll 



[2(5-5rcs)/r'J-'] +1 



(15) 



where a(5) — a{B) — i/3(5) and similarly for Orcs and 
flbg. The peak profiles for the elastic and total inelastic 
cross sections are given by combining these with Eqs. [3 
and[6l 

Some useful properties of the scattering length and 
cross sections follow from simple geometrical consider- 
ations. The complex scattering length a{B) has a value 
flbg far from resonance and describes a circle of radius 
ai.os|/2 in the complex plane as B is tuned across the 
resonance. If we write the resonant scattering length as 



I exp (2i6'i.cs) : 



(16) 
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FIG. 4: The resonant circle in the complex scattering length, 
showing the extent to which the imaginary part of Ores can 
reduce inelastic scattering. 



where ^rcs = 0o + ^ockbg, then the circle is as shown in 
Fig. m The smallest value achieved by P{B) is 

/3min = Phg - ^|arcs|(l - COs26',.es), (17) 
which occurs -Smin ~ -^rcs ~t~ ^min 

r'^° /2 with 

1 /'^ 

Xmin = -(arcs/Acs) " [(arcs/Acs)^ + l] " • (18) 

This defines the smallest value of the total inelastic cross 
section through Eq. [5l Unitarity requires that 1 5*001 < 1 
and P{B) > 0, so the limits on the possible values of O^cs 
are 

cos 26,,, > cos = 1 - ^ . (19) 

l^rcs I 

An obvious special case of this is that, if Obg is real, 
0res = so that Ores IS also real. 

Some examples of the possible behavior are illustrated 
in Fig. [5] for a case with moderately strong background 
inelasticity, /3bg — |arcs|/2 and |abg| ~ 2|arcs|- For these 
parameters, cos6'™g'^ — 0. When 6i-cs is close to its maxi- 
mum value, P{B) dips close to zero and the total inelastic 
cross section shows a trough that can reduce inelastic col- 
lision rates by more than a factor of 10. The elastic cross 
section also oscillates, but if |aios| <C l^bgl the oscillation 
is relatively weak and if abg and P,cs have opposite signs 
it is peak-like rather than trough-like. 

II. RESONANCES IN HE + NH (n = 1) 

We have carried out bound-state and scattering cal- 
culations on He -I- NH (n = 1) in a magnetic field us- 
ing methods almost identical to those used previously 
for n = [1^. The bound-state Schrodinger equation 
was solved using the BOUND package [3^, as modi- 
fied to handle magnetic fields |29(|. In the presence of a 



magnetic field, the total angular momentum is no longer 
a good quantum number. The calculations are there- 
fore carried out in a completely decoupled basis set, 
\nmn)\sms)\LmL), where s = 1 is the electron spin of 
NH and L is the end-over-end rotational angular mo- 
mentum of He about NH. All the m quantum numbers 
represent space-fixed projections on the axis defined by 
the magnetic field. The only good quantum numbers are 
the parity (— l)"+^+i and the total projection quantum 
number Mtot — "m-n + "m-s + wl- 

BOUND propagates a set of coupled differential equa- 
tions outwards from a point i?min, deep inside the inner 
classically forbidden region, and inwards from a bound- 
ary point i?max at long range. The two solutions meet 
at a matching point i?mid, and bound-state eigenvalues 
are found by locating values of the energy for which the 
inward and outward solutions can be matched. The pro- 
cedure used by BOUND is to seek energies at which one 
of the eigenvalues of the log-derivative matching matrix 
is zero [42j. 

For true bound states, i?max can be placed in the outer 
classically forbidden region. However, in the present 
work we are dealing with states of He-NH that lie close 
to the n = 1 threshold and are thus more than 30 cm~^ 
above the n — thresholds. Since there are open chan- 
nels, these are actually quasibound states and can pre- 
dissociate to form He -I- NH (n = 0). Nevertheless, they 
can still be located by artificially applying a bound state 
boundary condition at i?max, and this is how BOUND is 
used in the present work. 

Applying a bound-state boundary condition has the 
side-effect of box-quantizing the continuum states above 
both the n = and n = 1 thresholds. The resulting 
artificial bound states are easily identified because their 
energies depend on i?max- Difficulties arise only if an 
artificial bound state lies accidentally close to the level 
of interest, in which case the two states can perturb one 
another. Fig.l^lshows an example of artificial levels cross- 
ing the real levels as a function of i?max- For the case of 
He-NH it was usually possible to estimate the positions 
of the physical bound levels to within 0.01 cm~^. The 
perturbations are a measure of the genuine couplings to 
the continuum and are comparable to the width of the 
quasibound state, so this accuracy is sufficient for use in 
locating resonance positions. 

Fig. [7] shows the quasibound states of He-NH near the 
n = 1 threshold with artificial levels removed. Crossings 
between quasibound states and thresholds will produce 
zero-energy Feshbach resonances in s-wave scattering if 
an L = scattering channel is permitted by the con- 
straints on parity and Mtot- For A/tot — this occurs 
only for thresholds corresponding to rrij = 0, as shown 
by the circles in Fig. [71 

Each NH monomer level with n = 1 is split into 3 
components with j = 0, 1 and 2 by coupling to the spin 
s = 1. The He-NH levels closest to the n = 1 thresholds 
have predominantly L — 2 character. At zero field the 
total angular momentum J is a good quantum number, so 
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FIG. 5: (Color online). Some examples of the possible behavior for a case with moderately strong background inelasticity, 
/3bg = |arcs|/2 and |abg! = 2|ai.os[. The left-hand side shows the real and imaginary parts of the scattering length and the 
right-hand side shows the elastic and total inelastic cross sections (with inelastic cross sections calculated for a wave vector 
ko = 10'^ \ (2rcs I )• The 4 panels show (from top to bottom): (i) positive ttbg; ^rcs — — 20 ; (ii) positive o^bg; ^res — 44°; (iii) 
negative abg, 6^108 = —44°; (iv) negative Obg, ^ros ~ —20°. 
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FIG. 6: The pattern of levels from bound-state calculations 
on He-NH near the n = 1 threshold as a function of iimax- 
The calculations are for B — 4000 G, even parity, Mtot = 0. 
Note the artificial bound states with energies that decrease 
with 7?max, crossing and interfering with the physical states 
at constant energy. The lower panel shows a more detailed 
scan across a small region of -Rmax showing avoided crossings 
between real and artificial states. 

for He-NH each {n, j, L) level splits into min(2j'+l, 2L+1) 
components with different values of J. These splittings 
are barely visible in Fig. [71 so Fig. [S] shows an expanded 
view of the levels corresponding to {n,j,L) = (1,1,2) 
for all allowed values of Mtot- It niay be seen that the 
zero-field levels with J = 1, 2 and 3 are split by about 
0.04 cm"^. When a magnetic field is applied, each level 
splits into 2 J -I- 1 components with different values of 
Mtot. The J quantum number remains a useful label for 
magnetic fields up to about 200 G, but above that the 
levels of different J are strongly mixed. By about 600 
G the levels have separated into 3 groups that may be 
labelled with an approximate quantum number mj that 
takes values -fl, and —1. The levels corresponding to 
(n, j, L) show similar but more complex behavior. 

Once the crossing points have been located in Fig. [3 
the next stage is to carry out scattering calculations, 
holding the kinetic energy fixed at a small value (10~^ 
K in the present work) while sweeping the magnetic 



FIG. 7: (Color online). The pattern of levels from bound-state 
calculations on He-NH near the n = 1 threshold, with arti- 
ficial levels removed, as a function of magnetic field B. The 
calculations are for even parity, Mtot = 0. The circles show 
crossings between bound states and thresholds that produce 
zero-energy Feshbach resonances. 
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FIG. 8: (Color online). The pattern of He-NH levels arising 
from (n,j, L) = (1,1,2) as a function of magnetic field B. 
The calculations are for even parity and all allowed values of 
Mtot. 
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field across the resonance. This was done using the 



MOLSCAT package [3l|, as modified to handle magnetic 
fields [29|. MOLSCAT solves the Schrodinger equation 
by propagating a set of coupled differential equations out- 
wards from i?min to a matching point i?max at long range 
using basis sets and algorithms very similar to BOUND. 
The major difference is that it carries out only outwards 
propagation and matches to scattering boundary condi- 
tions at i?max, so that there is no artificial quantization 
of the continuum. 

MOLSCAT produces the S matrix and eigenphase sum 
S for each magnetic field B and evaluates the correspond- 
ing scattering length from Eq. |4l The next step is to 
fit to Eq. [7]to obtain the resonance position B^^s 

and width r^. This is done with the RESFIT pack- 
age (43j , which includes a quadratic polynomial in B for 
the background term Ebg. RESFIT then proceeds to fit 
the individual diagonal S-matrix elements to Eq. [8l hold- 
ing i?rcs and Tb fixed at the values obtained from the 
eigenphase sum. It represents the magnitude and phase 
of each Su with a quadratic polynomial and produces a 
complex number g^^ that describes the resonant circle in 
each S-matrix element. This fit provides all the param- 
eters required to define the real and imaginary parts of 
Obg and Ores needed for Eqs. [TT] and 1151 In practice we 
use a constant background term (obtained by evaluating 
the background polynomials a,t B = -Bics) in plotting the 
results of the formulae below. 

The results of fitting parameters to several resonances 
are shown in Table [J The first point to notice is that 
the resonances are all very wide, with iTsl > 30 G. This 
contrasts with the n = resonances previously charac- 
terized i29i], which had iFsl < lO'^ G. The difference 
arises because the n = 1 closed channels involved here 
are directly coupled to n = channels by the (weak) po- 
tential anisotropy, whereas the n ~ closed channels in- 
volved in our previous work were only indirectly coupled 
to open channels by a second-order mechanism involv- 
ing both the potential anisotropy and spin-spin coupling. 
The second point of interest is that arcs has a significant 
imaginary part in each case, with /S^cs/ctres considerably 
greater than Phg/ctbg- Because of this, there is signif- 
icant asymmetry in the calculated resonant line shapes 
for a{B) and f3{B). 

The real and imaginary scattering lengths for the res- 
onance labelled 1 in Fig. [7| are shown in the upper panels 
of Fig. [9] In each case the black points shows the numer- 
ical results from MOLSCAT and the red line shows the 
result of Eq.[TT]with the assumption that Oi-es is real (i.e., 
Ores replaced by |aros|)- It may be seen that there are sig- 
nificant discrepancies, shown by the red lines in the two 
lower panels. By contrast, if Oics is allowed to be complex 
(Eq. [TS]) , we obtain almost perfect fits to the numerical 
results except for the neglect of the field-dependence of 
the background scattering length Obg. The result is too 
close to the points to show usefully in the upper panels of 
Fig. m but the difference between the MOLSCAT results 
and those given by Eq. are shown as the blue lines in 



Figs. 1 

These results verify that Eq. [15] gives a correct account 
of the the behavior of the scattering length in the pres- 
ence of significant (but still small) background inelastic 
scattering. The complex nature of aj-os manifests itself 
in a slight tilting of the circles in the S-matrix elements, 
as shown in Fig. [31 The tilt is too small to be shown 
graphically, but may be seen in the numerical values in 
Table m 

The asymmetries seen in Fig. [Sj are much smaller than 
the possible asymmetries shown in Fig. [5l but they nev- 
ertheless serve to illustrate the principle that resonant 
signatures can be asymmetric and that inelastic cross 
sections can show troughs as well as peaks near reso- 
nance. In future work we will investigate systems with 
stronger background inelasticity in which more dramatic 
asymmetries can be expected. 



A. Computational details 

The bound and scattering calculations in the present 
work used a basis set with Umax = 5 and imax = 5. 
The coupled equations were solved using Johnson's log- 
derivative algorithm [3] with i?min — 1.7 A and a 
step size of 0.025 A. The bound-state calculations used 
^max = 12 A except where stated otherwise, and the 
scattering calculations used i?max — 100 A. 



III. CONCLUSIONS 

We have investigated the behavior of low-energy Fesh- 
bach resonances for a case where there is significant back- 
ground inelastic scattering far from resonance. We have 
located low-energy Feshbach resonances in s-wave scat- 
tering of NH ('^S^, n = 1) with He. The resonances are 
much wider (F^ > 30 G) than those found previously 
for scattering of NH (n = 0) (Fs < lO'^ G) We 
have investigated how scattering lengths and cross sec- 
tions vary as a resonance is tuned across threshold using 
an external magnetic field. 

The strength of a resonance can be characterized by 
a resonant scattering length Ores- If only elastic scatter- 
ing is possible, ares is infinite and the scattering length 
passes through a pole as a bound state crosses thresh- 
old. We have shown previously [1^ that this behavior is 
modified in the presence of inelastic collisions and that in 
some cases the scattering length shows a small oscillation 
rather than a pole. 

The key result of the present paper is the demonstra- 
tion that Ores can be complex rather than real and that 
this allows both the real and imaginary parts of the scat- 
tering length (and thus elastic and inelastic cross sec- 
tions) to show both peaks and troughs near resonance. 
We have shown that the real and imaginary parts of the 
scatte ring length follow analytical formulas given previ- 
ously |33 |. 
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TABLE L Parameters of magnetically tuned Feshbach resonances in He + NH (n — 1) collisions at -Ekin = K, corresponding 
to fco = 3.2189 X lO"* A"^. Resonances 1 to 3 correspond to the correspondingly numbered circles in Fig. [T] while resonance 
4 is for a different value of Mtot- 
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FIG. 9: (Color online). Upper panels: the real (left) and imaginary (right) parts of the scattering length for resonance 1 as a 
function of magnetic field. The black dots show numerical results from MOLSCAT and the red line shows the result of using 
Eq. (|ll|l with a real value of arcs- Lower panel: deviations between the numerical results and the formula using real (red, Eq. 
111(1 and complex (blue, Eq. [15} values of arcs. 



The present work offers the hope that tuning close 
to a Feshbach resonance can be used to reduce inelastic 
collision rates and thereby allow evaporative or sympa- 
thetic cooling in cases where collisional trap losses would 
otherwise prevent it. In future work we will investigate 
whether additional influences such as electric fields can 
be used to control the maxima and minima and allow 
detailed control of collision rates. 
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